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^3 ■ We show that the popular ILC approach is unstable in respect to the division 

of the sample of map pixels to the set of "homogeneous" subsamples. For suitable 
choice of such subsamples we can obtain the restored CMB signal with amplitudes 



o 



ranged from zero to the amplitudes of the observed signal. We propose approach 



u: which allows us to obtain reasonable estimates of Ct at £ < 30 and similar to WMAP 

Cg for larger £. With this approach we reduce some anomalies of the WMAP results. 
In particular, our estimate of the quadrupole is well consistent to theoretical one, the 
effect of the "axis of evil" is suppressed and the symmetry of the north and south 
galactic hemispheres increases. This results can change estimates of quadrupole 
polarization and the redshift of reionization of the Universe. We propose also new 



simple approach which can improve WMAP estimates of high £ power spectrum. 



ON 

o 

00 

O " INTRODUCTION ing methods [Ulll^, harmonic ILC [l3fl, fast 

O ! 

independent component analysis (FASTICA) 



14j etc. Among these approaches the ILC 



During last years fundamental results are 

obtained with the analysis of fluctuations of method is ver y convenient because in fact & 

relic radiation |Jfl observed by WMAP mis- rec l uires minimal additional assumptions in 

sion. Key problem of such analysis is the res P ect to the se P arated si S nals - Detailed 

cosmic microwave background (CMB) com- di scussion of the ILC approach with many 

ponent separation from the Galactic fore- corrections can be found in U, llQfl . The in- 

grounds in the pixel domain. Several ap- stability of the low multipoles reconstruction 

preaches were used to separate CMB from with the ILC method owin S to the correla - 

the observed signal. They are internal lin- tion between l M CMB and foregrounds was 

ear combination (ILC) and maxima entropy discussed in M- Recently some problems 

methods Q, Q, the blind and Wiener filter- arisin S with the ILC method were discussed 
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16|. CMB map and C e . 



In Planck review [17[, there is considered 
the final component separation pipeline for 
the Planck mission, which involves a com- 
bination of methods and iterations between 
processing steps targeted at different objec- 
tives such as diffuse component separation, 
spectral estimation, and compact source ex- 
traction 

At the same time some anomalies in re- 
sults of WMAP team are widely discussed. 
Among other these are the small amplitude 
of quadrupole component, unexpected corre- 
lations between components with £ = 2 & 3 
("axis of evil"), noticeable asymmetry be- 
tween north and south galactic hemispheres, 
existence of few deep walls in the CMB map 
etc. Final step with these discussions is paper 

18[ where all these anomalies are explained In Section 2 we represent four different 
as random fluctuations. procedure which can be used for the divi- 

Special problem is the analysis inhomoge- sion of the map pixels in the set of 'homo- 
neous map for which the amplitudes of fore- geneous' subsamples with analytical and nu- 
grounds strongly vary over the map. In this merical estimates of efficiency CMB compo- 



In this paper we show that the ILC 
method is unstable in respect to the defini- 
tion of " homogeneous" regions. As is shown 
below different criteria of homogeneity and 
corresponding division of the full sample of 
map pixels to set of " homogeneous" subsam- 
ples leads to different CMB maps and even 
different Ce. Thus, for suitable procedure we 
can obtain the CMB signal in wide range of 
its amplitude. In fact these amplitudes can 
vary from zero to the amplitude of observed 
signal. 



case the analysis 
as is described in 



3ecomes more complex and nent separation. In Section 3 we apply our 
4J it includes the division of "best" approach to the observed Q and V 



the map to set compact more homogeneous channels of WMAP and show that we can 

regions for which the component separation suppress some of the anomalies noted above, 

is performed independently. However in 4j Sec. 4 includes the summary of our results 

the choice of 12 such regions is not uniquely and discussion of methodical problems. In 

determined. Different definitions of "homo- particular, we propose new approach for the 

geneity" of selected subsamples are possible analysis of high £ power spectrum which can 

what leads to different final estimates of the improve now available results. 
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SEPARATION OF THE CMB SIGNAL 
WITH ILC APPROACH 

The ILC approach 

The observed map is builded set of 
pixels each of which contains combination 
S(6i) of the CMB signal C(^) and the fore- 
ground F(8i). If we have maps at two differ- 
ent frequencies then we can write 

Si{6i) = C{6i) + Ftfi) , (1) 

s 2 (e i ) = c{e i ) + F 2 {e i ), 

and we like to perform the linear extraction 
of the CMB signal as follows 

c(e i ) = as 1 (e i ) + (i-a)s 2 (e i ) 

= S 2 (6 i ) + a[S 1 (6 i )-S 2 {e i )], (2) 

The general expression for a determined 
by the condition of minimal dispersion of 
cleaned map is 

a = -(Q 2 Q 12 )/(Q 2 12 ), (3) 

o* = (C 2 ) - (C) 2 = (Ql) - (Q 2 Q 12 ) 2 /(Q 2 12 ) . 
Here 

Q 1 {9 i ) = S 1 {e i )-{S 1 ), Q 2 (9 t ) = S 2 (9 t )-(S 2 ) 

QxiiOi) = QiW-Q^Oi), (Qi) = (Q 2 ) = . 

and () means the averaging over the consid- 
ered subsample of pixels. 



However, as is seen from (P&[2]), 

a((l - Fx/Fa)) = 1, a = a f = -(l- {F 1 /F 2 )\A) 

where in accordance with the main ideas of 
the approach we consider a as a constant. 

Relation (TJJ points out the best value of 
the parameter of separation a = otf. This 
value depends upon the ratio F\jF 2 and the 
scatter of a is determined by the scatter of 
this ratio for the subsample used. Moreover, 
two values of the parameter of separation, 
a (EJ) and a/ (jlj), are different and this dif- 
ference decreases for decreased scatter of ra- 
tio Fi/F 2 . This means that in order to im- 
prove the separation we must divide the full 
sample of pixels to set of more homogeneous 
subsamples using the distribution of ratios 
F\0j) I F 2 (8i) . After the component separa- 
tion within these subsamples we get set of 
cleaned pixels sum of which forms the cleaned 
map and allows to perform further analysis of 
this map with better precision. Example of 
such component separation is considered be- 
low (model 1). 

However, such approach cannot be used 
in practice when the foregrounds are a pri- 
ory unknown and for the component separa- 
tion we would have to use criteria expressed 
through the observed signals. As we show be- 
low the cleaned map strongly depends upon 
these criteria. 

In the further analysis we consider the pix- 
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els as independent ones and ignore the pos- In main this approach is similar to that 
sible correlations of the signal amplitude in used in [4J in order to take into account the 
the neighboring pixels. The inclusion of such inhomogeneities of the foreground. However, 
correlations allows to improve the component their selection of 12 pixel subsamples differs 
separation but makes the procedure of sepa- from ones discussed below. Our analysis con- 
ration more complex. firms that the correct result can be obtained 
As demonstration of these statements we only for the known a priory foregrounds. In 
consider below both analytically and numer- all other cases we can obtain the approximate 
ically four models of map division on "homo- estimate of the CMB signal only. But devi- 
geneous" subsamples prepared with various ations between the input and restored CMB 
definitions of " homogeneity" . We determine signals depend upon the criteria homogene- 
the "homogeneous" subsamples in respect to ity and decreases for less A. For larger A all 
the function G of amplitudes of signals approaches give comparable results. 

G i = G{6 i ) = G(S 1 {6 i ),S 2 (e i )) 

The i th bin contains K, pixels for which we Four models of separation of the CMB 

have 



signal 



i < d/A < i + 1 (5) 



The theoretical consideration reveals the 
where A is a given common width of the bins. main in fl uences Q f the selection criteria but 
The bin center is the mean amplitude of the real estimates of quality of separation can 
function G^ ^ e f ouno i with, simulations only. To test 



the various methods of component separation 



(G i )=Y,Gik/K i , (6) 

fc=i we generate the CMB signals with the stan- 

By the way for all bins we have the symmetric dard power spectrum and Gaussian distri- 

distribution of functions G lk with bution of amplitudes, using the foregrounds 



from WMAP 191 we transform the generated 
\S k \ = \G ik -(G i )\<A, (6 i ) = (G i -(G i )) = 0. ^ 

CMB signals to observed ones and separate 

For each subsample we obtain on accord- the CMB signals with various approaches. 

ing to the standard relation ([3]) and get the The final estimates of precision achieved for 

CMB signal, C(8ik) for each pixel of consid- the full map relate to the comparison of in- 

ered subsample with relation |2j troduced and restored Ct. 
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model 1 

Let us consider the set of subsamples with 

0(9,) = F 1 (9 i )/F 2 (9 i ) = 1 + P + 5(9,) , (7) 

(G) = l + P, \5\<A, 

F 1 = F 2 (l+P + 8 i ), (Q 12 ) = (3(F 2 ) + (F 2 5) , 

Here 1 + (3 is the center of the subsample and 
$i = S(9i) characterizes the (small) random 
scatter of the pixel amplitude in respect of 
the central point ((5) = 0) . 
For such subsample we get 



a 



1 + o(5) _1 

'P + o(sy af ~ 



and for 5 — >■ we have a — > aif — — 1/(3, 

C(9 i )=C(e i ) + A c (9 i ), 

A (a\ r?/a , (F 2 5)/(F 2 ) + 5(9 l ) 
Ac(ft) = F 2 (9,) oc5. 

As is seen from this relation 



(8) 



A c (0i) ^ for 5^0 



(9) 



For such choice of the pixel subsamples we 
get accurate component separation precision 
of which depends upon the bin size, A, and 
increases for smaller A. Numerical simula- 
tions confirm this conclusion. 



where again 1+/3 is the center of the subsam- 
ple and 8(9i) characterizes the (small) ran- 
dom scatter of the pixel amplitude in respect 
of the central point ((5) = 0, \5\ < A). In 
the case 

SMS^) = s 2 (0O(0+<Si), Q 2 ( di ) = s 2 (ei)-{s 2 ) , 

Qi2(9i) = (3Q 2 (9 t ) + S 2 (9 t )5(9 t ) - (S 2 $)ll) 

(Q 2 12 ) = (3 2 (Q 2 2 ) + 2(3(5S 2 Q 2 ) + o(5 2 ) , 
(Q2Q12) = (3(Q 2 2 )+(5S 2 Q 2 ), a » -l/(l+o(5) . 
Therefore, 



cr. 



c 



(Ql) 



mi) + o(5) 

l + o 1 (5) 



oc o(5) . 



l + o 2 (5)_ 

Thus, we see that C(9/) oc A, a 2 c oc A, and 
for A ^ we have C(9i) ->■ 0, a% ^ 0. For 
such pixel subsamples we get the extremal re- 
sult - the signal CMB equal zero. The same 
result can be obtained for an arbitrary func- 
tion G = G(Si/S 2 ). Numerical models con- 
firm this tendencies. 

model 3 



model 2 

Let us consider the set of the pixel sub- 
samples with 



Let us consider the set of pixel subsamples 
with 

G l = S 1 (9 l ) = S [l + 5(9 l )], (13) 



Gi = S l (9 i )/S 2 (9 l ) = l + (3 + 5(9,) (10) 



Qi2 = S 5(9 t )-Q 2 . 



6 



Here So is the center of the subsample and in Q and V bands are presented in FigJTJ It 

5(9i) characterizes the (small) random scat- is interesting that the best reconstruction is 

ter of the pixel amplitude in respect of the obtained for the larger A and for restored and 

central point ((5) = 0, \5\ < A/So). In the input signals the ratio Ce/C in decreases with 



case 

(Ql 2 ) = (Ql)-2S (5S 2 ) + S 2 (5 2 ), 

(Q2Qi2) = -{Ql)+So(6Q 2 ), a= 1 ~° l{6) 

1 - o 2 (()J 

C{9 i ) = S 2 (6i)(l -a) + aS [l + 5(001(14) 

Thus, for 5 — > we get 

a ->■ 1 G[6i) ^ S , o%->0. (15) 

For such choice of the function Gi ( fl~3l) we 
get unexpected result - for small A — > the 
signal CMB is equal to Si = So. Numerical 
simulations confirm these tendencies and as is 
seen from the Table 1 for small A the selected 
signal C is quite close to the input one Si and 
strongly differs from S 2 . For larger A this 
difference disappears. 



A. 
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FIG. 1. For the model 3 the reconstruction 
of modeling CMB signal for the foregrounds 
in Q and V bands. Results are normalized 
on the input signal and are shown for A = 
2,0.2&0.002mif (solid, dashed and dot-dashed 
lines) . 



model 4 



TABLE I. Two examples of the reconstruction Let us consider the set of pixel subsamples 
of the CMB signal with the model 3 (arbitrary with 
units) 



-Gi = Sd 



A 



Npixels (S 



(S 2 



(C) 



Q 



12 



)-S 3 (9 i ) = F 1 (9 i )-F 2 (9i) = P(lJl 
(35(9,) (6) = 0, \6i\<A/l3. 



0.2mK 256129 8.7 ± 5.4 5.6 ± 5.3 3.3 ± 6.3 
0.002mK 2835 1.0 ± 0.57 -20. ± 17. 0.4 ± 0.8Here (3 is the center of the subsample and 
0.002mK 2922 3.0 ± 0.6 -18. ±18. 1.3 ±1.6$ = K^i) characterizes the (small) random 



scatter of the pixel amplitude in respect of 
For A = 2, 0.2 & 0. 002m.fr reconstruction the central point. As in the model 1, so de- 
of the modeling CMB signal with foregrounds termined function Gi depends upon the fore- 
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grounds only what is some advantage of this 
approach. In the case 

(Q 2 Q 12 ) = (3(S 2 5), (Qj 2 ) = (3 2 (5 2 ), a = - 

C{Bi) = S 2 {9 i ) - [1 + 5(^)](5 2 5)/(5 2 ) (17) 

4 = (Ql) - (S 2 5f/(5 2 ) 

For such choice of the function Gi results 
depend upon the bin size but even for 5 — > 
they are not tend to real CMB signal. In the 
case the choice of optimal A can be done with 
simulations. 

Examples of such reconstruction of the in- 
put CMB signal with A = 2, 0.2, & 0.002mA" 
are presented in Fig. [2j As is seen from this 
figure reconstructed signal is weakly sensitive 
to used small A and is oscillated around the 
level C e /C in ~ 1. - 1.1. 

The difference between models 3 and 4 is 
illustrated by Fig. [3] where we see the proba- 
bility distribution function P(a) for fraction 
of pixels versus the separation coefficient a. 

POWER SPECTRUM FROM Q AND V 
BANDS OF WMAP 

As was found in previous Section the best 
reconstruction of Ci is possible with approach 
used in the model 4. Applying this approach 
with A < 0.2mK for Q and V bands of the 
WMAP maps we get which significantly 
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FIG. 2. For the model 4 the reconstruction 
of modeling CMB signal for the foregrounds in 
Q and V bands. Results are normalized on 
the amplitude of input signal and shown for 
A = 2,0.2&0.002miT (solid, dashed and dot - 
dashed lines). 
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FIG. 3. For A = 0.02mK the distribution of the 
fraction of pixels vs. a for model 3 (top panel) 
and model 4 (bottom panel) 

differ from ones presented in WMAP publi- 
cations. In these cases we have from several 
tens to several thousands of 'homogeneous' 
regions instead of 12 regions used in WMAP 
analysis. These AT/ are plotted in Fig. 0] 
and a 2m are listed in Table [Til However, for 
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broad bins with A > lOmK our results be- In contrast, for our parameters of 
come quite similar to the WMAP ones. quadrupole we get 



TABLE II. Amplitudes of quadrupole compo- 
nents in [i K for A = 2fiK 

WMAP Model 4 



«2,0 


11.48 


-65.2 


(32,1 


-0.05 


-13.8 


0-2 -1 


4.86 


9.0 


^2,2 


-14.41 


-17.3 


(32,-2 


-18.80 


-11.0 



Ai= 68.3//AT, (/,&) = (-75°, 9.1°), 

A 2 = U.QnK, (l,b) = (13.1°, -8.(29) 

A 3 = -80A/iK, (I, b) = (60.°, 77.4°) , 
with 

AT 2 = -A( AiA2+AiA 3 +A2 a 3 ) = 1070/i^ 2 . 

The orientations ( |22l) differ from both the 
With listed in Table [III we get for the ^ 
quadrupole ((, 6) D = (-96°, 48°) , 

AT^IOTOM- 2 , (18) and from or . entations m 

what is close to theoretical expectations 

AT 2 ^ 1250/xK 2 . (19) 
and exceeds estimate obtained by WMAP 



team 



4] 



ATq « 249/iiT 2 . 



(20) 




lOO 



As is well known, the five quadrupole co- FIG. 4. The 10~ 3 AT 2 for the WMAP data 

efficients are equivalent to the components of (points) and obtained according to the method 

a symmetric traceless tensor. For the princi- used for the model 4 ( stars) with A = 2 o^K. 

pie values and orientation of tensor axes for SoM and daghed lines ghow the theoretical i y 



ex- 



the 3 years WMAP quadrupole we have |20| pected yalueg and their gcatter _ 
Ai = 27.1//AT, (I, b) = (-0.8°±13°, 63.3°±1°) 



As is seen from Fig. H] the most serious 

differences are found for £ = 2, 4, and for 
20)) 

even £ < 30. For these even £ our estimates 
A 3 = -40 [iK, (l,b) = (-77.6°±5°, 6.5°±4°pf C t exceed ones obtained by WMAP by a 



A 



12.9/xFT, (I, b) = (15.5° ±3°, 25.8° ±1 
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factor of ~ 1.5 what emphasizes the symme- 
try of the CMB signal in north and south 
hemispheres. For £ > 30 the difference be- 
comes small. It is interesting that for odd £ 
deviations from WMAP results regularly do 
not exceed 10%. The random scatter of the 
method depends upon the bin size used, A, 
but not exceed ~ 10 — 15% what does not 
distort essentially our estimates of the power 
spectrum. 

These results noticeably change low £ part 
of the power spectrum and significantly sup- 
press the effect of "axis of evil". However, 
they do not distort strongly main conclusions 
of WMAP which are weakly depend upon 
this part of the power spectrum. 

Let us emphasis only that new estimate 
of Gi can noticeably changes the estimates of 
the quadrupole polarization and, therefore, 
the redshift of reionization. 

SUMMARY AND DISCUSSION 

In this paper we show that the separa- 
tion of foregrounds and the CMB signal with 
the ILC method strongly depends upon the 
choice of 'homogeneous' subsamples of pixels. 
For foregrounds presented in WMAP papers 
our more stable estimates of the CMB fluctu- 
ations are obtained for the selection criteria 
used in the model 4. Theoretical consider- 
ation (Tl6]) shows that with this approach we 



cannot perform the very high precision clean- 
ing. However, numerical analysis demon- 
strates that for suitable choice of the bin size, 
A, the precision a ~ 10% can be achieved. 
It can be expected that the application of re- 
fined technique developed by WMAP team 
will allow to decrease the errors up to values 
presented in ^4 1 . 

Main results 

The best results are obtained for the fre- 
quency channels Q & V and are presented in 
Fig. HJ Main results of our analysis can be 
summarized as follow: 

1. The measured amplitude of quarupole 
is more than that given by WMAP 
by a factor of 2.1 what eliminates 
disagreement between the theoretically 
expected and measured values. 

2. The coordinates of the quadrupole are 
changed while our estimates of the oc- 
tupole remain the same as in WMAP. 
This fact substantially reduces the ef- 
fect of " axes of evil" . 

3. All even Ct with 4 < £ < 20 are more 
then those given by WMAP by a factor 
of ~ 1.5 — 2 what emphasizes the sym- 
metry of the CMB signal in north and 
south hemispheres. 
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4. Deviations of odd Ci from that given Methodical comments 
by WMAP do not exceed a factor of 

^2 The considered models allow us to obtain 

some inferences related to the method of lin- 
ear component separation. Thus, we see that: 

5. At i > 30 deviations of our estimates 

from the WMAP data do net exceed 1- The method of linear component sep- 
~ 5%. aration is unstable and the resulting 

CMB map strongly depends upon cri- 
teria homogeneity used for the selection 

6. At i < 20 the expected error of mea- 

of the set of subsample under consider- 

sured Cg is ~ 10%. 

ation. 



These results indicate that the main con- 
clusions of the WMAP team remain correct. 
However, the change of the large scale char- 
acteristics leads to the moderate change of 
estimates of <j$ and especially the estimates 
of low £ polarization and, therefore, the red- 
shift of reionization of the Universe. These 
corrections could be important for analysis 
of the epoch of reionization and formation of 
earlier galaxies. 

The further more detailed analysis of pos- 
sible divisions of the full sample of pixels to 
the 'homogeneous' subsamples can find more 
effective methods of subsample selection than 
that used in the paper. In particular, the ac- 
count of correlation of the signal amplitude 
in neighboring pixels can improve the quality 
of the cleaned map of the CMB signal. 

Of course, this approach can be extended 
for the three and more frequency channels. 



2. The best separation is possible with 
using the foreground measurements 
(model 1). However, such approach is 
of no concern for a practice as we do 
not know a priori the foregrounds. 

3. Models 2 and 3 demonstrate that with 
a suitable choice of the selection criteria 
we can obtain arbitrary estimates for 
the CMB signal. 

4. Reasonable estimates of the CMB sig- 
nal can be obtained with the selection 
criteria used in the model 4. How- 
ever even in this model the CMB signal 
can be found with errors which depend 
upon the bin size A ([5]) used for the 
subsample selection. 

5. Comparison of theoretical estimates of 
ac for models 3 and 4 with numerical 
estimates of Ci shows that sometime 
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the former ones do not characterize ad- 
equately the final precision achieved. 

It can be expected that final results de- 
pend upon the actual foreground. This infer- 
ence is confirmed by comparison results ob- 
tained for various pairs of frequency channels. 

Let us note that further cleaning can be 
performed by recurrent comparison of the 
cleaned maps obtained for two pairs of fre- 
quencies. With the WMAP data we cannot 
test this approach as the quality of maps ob- 
tained for QV channels significantly exceed 
the quality of maps found for other pairs of 
frequency channels. However, for many chan- 
nels of the PLANCK mission such approach 
becomes useful. 

Estimates of the high £ power spectrum 

As is well known for the real maps of the 
CMB with the finite number of pixels the de- 
termination of the power spectrum for larger 
I is complex because the polar regions with 
relatively small number of pixels along the az- 
imuthal coordinate cannot be used. By the 
way at high i we would have to analyze the 
noisy regions in the vicinity of equator what 
decreases the precision achieved. 

To decrease the influence of the noisy 
galactic equator we can use the simple proce- 
dure what is change of the map orientation. 
Indeed, if we will build the map in coordi- 









II 







FIG. 5. The CMB map for the Q channel after 
rotation of the coordinate system. 

nate system with the galactic equator situ- 
ated along some map meridian then we will 
have less noisy pixels situated along the map 
equator while some of the noisy pixels will 
be shifted to polar regions. Example of such 
map is presented in Fig. [5J 

Of course, such approach requires prepa- 
ration of two different maps one of which have 
the ordinary orientation and is used for the 
analysis of the low i part of power spectrum 
while second one with the orthogonal orien- 
tation can be used for analysis of high I com- 
ponents of the power spectrum. 

This approach seems to be quite effective 
but it must be tested with real repixelized 
maps. 
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